Quantum Tunneling and Scattering of a Composite Object: 
Revisited and Reassessed 
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This work presents an extensive exploration of scattering and tunneling involving composite ob- 
jects with intrinsic degrees of freedom. We aim at exact solutions to such scattering problems. 
Along this path we demonstrate solution to model Hamiltonians, and develop different techniques 
for addressing these complex reaction-physics problems, discuss their applicability, and investigate 
the relevant convergence issues. As examples, we study the scattering of a two-constituent deuteron- 
like systems either with an infinite set of intrinsic bound states or with a continuum of states that 
allows for breakup. We show that the internal degrees of freedom of the projectile and its virtual 
excitation in the course of reactions play an important role in shaping the S-matrix and related 
observables, giving rise to enhanced or reduced tunneling in various situations. 
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I. INTRODUCTION 



Reaction physics involving composite objects is a ma- 
jor and critical subject while encountered in the context 
of processes like fusion, fission, particle decay, as well as 
specific branches of science including chemistry, atomic 
physics, condensed matter physics, etc. In all these phe- 
nomena, more often than not, the scattering or tunneling 
object has its own degrees of freedom. Various pertinent 
scenarios have been explored earlier, where the tunnel- 
ing has been shown to be enhanced by the additional 
degree(s) of freedom, which may have arisen from the 
compositeness of the object [l|-[3|, from its interaction 
with another particle Q, or directly from quantum field 
excitations [9(. 

This is a complicated and generally non-perturbative 
problem, involving vastly different scales. While there 
are many techniques and methods of dealing with this 
problem, most of them involve simplifications. For ex- 
ample, some studies of the models that are similar to 
ours involve restriction on the range of energy of the pro- 
jectile [II], the mass-ratio of the constituents 0, the 
number of states available in the intrinsic system [3, Q , 
etc. In addition, most models exclude the possibility of 
virtual excitations of the object undergoing a reaction 
[US 0|- While simplifications work well at times, it is 
also common that the "slightly simplified" problem turns 
out to be very different from the original one. Moreover, 
some formally exact techniques, as demonstrated in Ref. 
[Tfl | and further discussed in this paper, do not necessarily 
provide a path to a convergent solution for an arbitrary 
subset of parameters. The paramount goal of this work 
is to find an exact solution to a given reaction problem, 
which is free from the above mentioned limitations and 
is reliably convergent. 



In order to reach this goal we limit our studies specif- 
ically to a problem in one dimension and to a composite 
object with two constituents, only one of which inter- 
acts with an external potential. A deuteron hitting a 
Coulomb barrier could be a fair example of such a pro- 
jectile. This picture has been modeled in several different 
ways in our work. However, the techniques that we de- 
velop and the study of how they work are general and 
are not limited, in their applicability, to our examples 
or models only. Moreover, we believe that many of our 
findings are generic and there are realistic situations that 
can be represented by even these simple models (ll) . 
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Our discussion is organized as follows. In Sec. [TTI we 
start by identifying our models and invoking some defi- 
nitions of reaction physics. Then in Sec. IIIII we examine 
a particularly simple example of a deuteron-like system 
reflecting from an infinite wall. This case provides an ex- 
cellent illustration of the pivotal role of virtual excitations 
in the dynamics. It also shows how the formally exact 
method of projecting the reaction dynamics onto the in- 
trinsic shell-model-like space could fail to yield reliable 
results. We put forward and demonstrate the Variable 
Phase Method (VPM) in Sec. [Tvl followed by solutions 
to various examples in Sec. [V] While the VPM has 
been used by others in the past, we extend it so as to 
include virtual channels. This novel extension requires 
us to explore the role of virtual channels and to discuss 
the convergence of solutions with the number of virtual 
excitations included in consideration. This is done in 
Sec. IVI1 A study of scattering and breakup of a system 
with a continuum of states is presented in Sec. IVIII The 
summary and conclusions are laid out in Sec. IVIIII 
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II. GENERAL DESCRIPTION OF THE 
PROBLEM 

A. The model 

Throughout this text we examine a one-dimensional 
problem. We consider a projectile which is a composite 
object made up of two particles which have masses mi 
and ?7i2 and are bound by an intrinsic potential v(xi — x%), 
where the particle coordinates are x\ and X2, respectively. 
This composite system interacts with an external poten- 
tial V{x\,X2)- The usual center-of-mass and relative co- 
ordinates are 



X 



miXl + m 2 X2 

M 



X = X\ — X2, 



(1) 



and the corresponding total and reduced masses are 

TOi7B2 



M = mi + m 2 , 



m 



mi + m 2 

The Hamiltonian for the system can be written as 

h 2 d 2 



(2) 



H 



2M dX 2 
where the intrinsic Hamiltonian 

h 2 d 2 



V(x 1 ,x 2 ) 



2m dx 2 ^ ^ 



(3) 



(4) 



has eigenstates ip n (x) with the corresponding energy 
eigenvalues e n : 



hip n (x) = e n tp n (x), n = 0,1,2,.. 



(5) 



The general scattering problem can be formulated in 
a traditional way, using the asymptotic form of the wave 
function. At X — > =Foo it is given by 



e lK *' x ^ R 



d>(A» 



\fWn 
nEopen V i '* 



E T 



nGopcn VT^nl 



nri iK n X 



(6a) 

<P n (x) (6b) 



respectively. The wave function above corresponds to an 
incident beam coming from the left with a particle in the 
intrinsic state (channel) n'. Here K n is the center-of-mass 
momentum of the system with total energy Et while in 
channel n, 



1 



K n (E T ) = -y/2M{E T -a r 



(7) 



Here, and later in this paper, the symbol '~' is used to 
indicate an asymptotic equality that involves only open 
channels n with Et > £ n - Contributions from the closed 
channels, for which Et < £ n , decay exponentially with 
distance from the scattering potential and are not present 



in the asymptotic form. The coefficients R and T are 
referred to as reflection and transmission amplitudes due 
to their physical meanings. |i? ran '| 2 and |r nrt '| 2 represent 
the probabilities for the incoming beam in channel n' 
to reflect and transmit, respectively, in channel n. The 
conservation of probability hence implies 



E 

nGopcn 



(\Rnr 



\T nn A z =1. 



(8) 



It should be mentioned that for the scattering problem 
to be fully determined one should consider, in addition 
to in]), an incident beam coming from the right, which 
gives rise to another set of reflection and transmission 
amplitudes. If and when it is necessary to distinguish 
between these two, the amplitudes in Eqs. ^ for the 
incident beam coming from the left and traveling in the 
positive x direction are denoted by R + and T + instead 
of just R and T; i?_ and T_ are used when they are 
associated with an incident beam traveling from the right 
to the left. 



B. The S-matrix 

While it is convenient to use the reflection and trans- 
mission amplitudes, the formal S'-matrix is still essen- 
tial for establishing a relation between this description 
and the traditional scattering theory. In addition, S'- 
matrix allows one to utilize the symmetries of the prob- 
lem, and determine relations among the amplitudes. De- 
spite S'-matrix being a textbook subject (l2l . [l3| . there 
are a few non-trivial features that emerge in the case of 
coupled-channel problems and with non-symmetric po- 
tentials |T^ - [l6| . We review some of them in what follows. 

Let us first consider the case of one open channel, 
which is of particular importance for many examples 
considered in this work. For a real potential barrier 
the transmission amplitude is symmetric between the in- 
coming beam traveling from the left and that from the 
right, which follows directly from the complex-conjugated 
Schrodinger equation, showing time-reversal invariance. 
The S-matrix can be defined in several different ways 
[14| . It is quite common to select a basis with incoming 
and outgoing waves so that S = 1 at high energies or 
in the absence of a potential barrier. An alternative ap- 
proach is to choose the S-matrix to be symmetric, which 
is possible because of the time-reversal invariance. Un- 
fortunately, it is impossible to accommodate both prop- 
erties simultaneously. We choose the second alternative 
and define the S-matrix using the following symmetric 
and asymmetric (in space) asymptotic forms of the in- 
coming waves at |A| — > 00, 

$+(A) ~ -^=exp(-iK\X\) and 
$-(x)^-^^-ex P HAlA|). 
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The outgoing-wave basis comprises the corresponding 
complex-conjugated wave functions. Then the S-matrix 
in terms of reflection and transmission amplitudes is 

q _ _l f (R+ + R-) + 2T (R--R+) \ 
2^ (R--R+) (R+ + R-)-2T J ■ 

Note that the S-matrix is symmetric, and T + = T_ = T 
due to time-reversal invariance. From unitarity of the 
S-matrix we find that = \R + \, \R%\ + \T%_\ = 1, and 
5ft [T*(R + + =0. The convenience with the above 
definition is that, for a symmetric potential, i?_ = 
Also, parity is a good quantum number, and hence the 
S'-matrix is diagonal with matrix elements 

S ± = -(R±T)= exp(2i5 ± ) . (9) 

An extension of the above discussion to a more gen- 
eral multichannel case is straight- forward [HI, [l?} • From 
unitarity it follows that 

R±R± + t±t ± = i, r\t t + tIr t = 0. (10) 

Time-reversal invariance leads to 

T± = T T , R± = R t ± . 

Finally, reflection symmetry of the scattering potential 
leads to R + = VR-V and T + = VTJP. The equalities 
for one channel are modified due to the different parities 
of the intrinsic states of the composite object. V denotes 
the parity operator in the channel space, so that V 2 = 1 
with Vnn' = S nn 'TT n , where ir n is the parity of the intrinsic 
state tp n (x). 

In the presence of reflection symmetry it is sufficient 
to consider only beams originating from the left and thus 
to deal only with R + and T + . From this point onward 
we omit the subscript + (hence returning to the nota- 
tions used in Eqs. ([5]). The symmetries discussed for the 
multichannel case are summarized as follows: 

R = R T , T = TT T T 1 (11) 



III. THE PROJECTION METHOD: EXAMPLES 
AND LIMITATIONS 

Before we actually describe the Projection Method, let 
us emphasize one important issue. While the observed 
picture (the S'-matrix for example) is seen through the 
asymptotic forms of the wave functions in the open chan- 
nels, the crucial dynamics of a scattering process takes 
place in the vicinity of the scatterer, and involves virtual 
(or closed) channels just as much as the open channels. 
The virtual channels are populated in accordance with 
the time-energy uncertainty, and lead to an immensely 
complicated process. Excluding the virtual channels from 
consideration could therefore lead to erroneous results in 
the observed quantities. A "simple" model example dis- 
cussed below elucidates both the complexity and the im- 
portance of virtual excitations. This model constitutes 
an infinite wall as a scatterer which interacts with only 
one of the two constituents of the composite object. We 
refer to this model as the "deuteron and Coulomb- wall" 
model, as defined in Sec. MI Al 

It is noteworthy that numerous authors [l8l - l20j . have 
worked on this subject, but reports of the findings are 
scarce. Mathematical issues, difficulties with stability of 
the solutions, and lack of appreciation from the scientific 
audience are some of the possible reasons. 

In the following method, referred to as the Projection 
Method, a solution is attempted by projecting the reac- 
tion dynamics onto the intrinsic basis set. In some way 
this approach is similar to various projection techniques 
used in nuclear many-body studies that involve reactions 

pis. 

We start our presentation by returning to the 
"deuteron and Coulomb- wall" model and to the Projec- 
tion Method. We draw some conclusions regarding the 
earlier discussions (l8l - [20l [23 - [27| by presenting an exact 
solution, showing limitations of the Projection Method, 
and highlighting the overall importance of this example 
for the understanding of reaction dynamics and develop- 
ment of techniques. 



RtR + TiT = l, R\VT) + {PT) jf R = 0. (12) 
We define 

S ± = ~{R±VT) (13) 

in this case, so that the S-matrix is symmetric, and the 
phase shifts approach zero in the limit of zero energy, 
since R = — 1 and T = in this limit. Conditions at 
other thresholds are related to Levinson's theorem which, 
for one-dimensional scattering, is discussed in Ref. (l5j . 



A. The deuteron and Coulomb-wall model 

For all the models in this work we assume that in the 
composite projectile, loosely referred to as the deuteron, 
only the second particle interacts with the potential, 
V(xi,X2) — > V{x%). In this section we concentrate on 
an example where the potential is represented by an in- 
finite wall or, the "wall": 

V(X ) — I °° W ^ ei1 < ^2 

I otherwise 
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The traditional textbook methods prescribe looking for 
a full wave function in the form 



where matrix D is defined as 



$(X,x) 



E 



Rn 



-iK n X 



(14) 

with the boundary condition x) — at X2 = 0. In 
contrast to the asymptotic form in ([6]), where the sum- 
mation includes open channels only, the sum here is over 
all channels and the expression holds for all values of 
X2 < 0. The meaning of the reflection amplitudes R nn i 
is therefore extended to include the virtual channels as 
well as open. The asymptotic form © is recovered at 
\X\ — > oo because the term corresponding to each vir- 
tual channel, say n, decays with distance (from the wall) 
through the exponential factor e~'^" x ', and therefore 
does not appear in the asymptotic sum. This exponent 
can be expressed in a generic way as e lKn \ x \ by assuming 
the principal branch of the square root in Eq. ([7} . The 
branch being specified allows one to consider momentum 
in a complex plane. 

The location X2 = in the boundary condition trans- 
lates into x = X\ and X = fax, as follow from Eqs. fTJ). 
Here we define relative masses as 

Mi,2 = mip/M, and fx = m/M. (15) 

Thus, the equation to be solved is 

^(mx, x)=0 (16) 

for all x's. 

The length scale for this problem is determined by 
a quantity A that is associated with the characteristic 
width of the intrinsic potential v(x). The intrinsic sys- 
tem also defines the energy scale, based on the usual 
coordinate-momentum uncertainty, as 



(17) 



2mA 2 ' 



In what follows we use A and e as our units of length and 
energy respectively. This is equivalent to using dimen- 
sionless energy units rescaled to e, namely, e n — > e n /e 
for the intrinsic energies, and E — > E/e for the center-of- 
mass kinetic energy ; and lengths rescaled to A, namely, 
x — > x/X. X — >• X/X and K n — > K n X for coordinate and 
momentum variables. Thus, it is assumed that A = 1 and 
e = 1 unless otherwise stated. The center-of-mass kinetic 
energy for an incident beam in channel 7i is E — Ef — £ n ^ 
in almost all our examples the incident beam is in the 
ground state channel, and therefore n = 0. 

Truncating the number of channels at some large N 
and looking for a solution in the space spanned by the 
functions ip n for n < N constitutes the projection ap- 
proach. Thence emerge the equations 



E 



Din [-im(K n > + K n )] 

\fWn~\ 



Rn 



(18) 



Di n (x) = / ip*(x) e* x ip n (x) dx. 



(19) 



This is the expectation value of the momentum shift op- 
erator in the intrinsic basis. Eq. (|18p is obtained by 
projecting the boundary condition onto the intrinsic ba- 
sis set. Note that for a virtual channel, the argument of 
the Z?-matrix becomes real and positive. 

At beam energies below the first threshold, when only 
the ground state channel (n = 0) is open, Eq. (fTl))) is 
particularly simple, since scattering is characterized by 
only a single phase of the reflection amplitude. Due to 
unitarity |i?oo| — 1; and the single S- matrix phase S is 
defined through e 2lS — — i?oo- Equation (JTHJ) then reads 



ijjo(x) sin [fiiKQX — 5] 



E 

nG closed 



R, 



nO 1 



\Kn\ 



where R' n0 = R n oe 



-»(iJ+7r/2) 



is real for any n. 



(20) 



To further illustrate the situation let us review two 
specific examples of intrinsic potential v(x) where the 
eigenstates ([5]) and the shift matrices (p~9|) can be found 
analytically. 



1. Infinite square well ("well") confinement 

In the first example, which is that of an infinitely-deep 
square well intrinsic confinement, the length scale A is 
defined so as to set the width of the well to 7rA, thus 



v(x) 



when |x| < nX/2 
oo otherwise 



(21) 



The eigenstates and the corresponding energies for this 
square well are 







7T\ 1 


sin 


(x -\ 


— n 






2/ J 



Sn-i = n 2 , (22) 



where n — 1, 2, . . . , so that the indices for both tp and 
£ start from 0, and the energy scale p7|) is e = £q. The 
corresponding shift matrix is 



Ann' ' k 



D, 



(-!)"+"' exp(^) 



exp I 



2 ) 



7T \(n + n') 2 + x 2 ] [(n - n 1 ) 2 + x 2 ] 



(23) 



2. Harmonic oscillator ("HO") confinement 



One could criticize the infinite square well potential as 
being too sharp and therefore leading to nonphysically 
high intrinsic excitations. Therefore, the harmonic oscil- 
lator intrinsic confinement v(x) — mw 2 i 2 /2 is presented 
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as a second example, which does not have this controver- 
sial feature. 

The unit of length h ere is d efined by the standard os- 
cillator length, A = y/ 6/raw. The eigenstates and the 
eigenvalues are defined in terms of the usual Hermite 
polynomials H n , 

ip n (x) = —j===H n (ir)exp ( ) , e n = (2n + l), 



where n = 0, 1, . . . , and the energy unit is e 
The corresponding shift matrix is 



(24) 



Am' (x) 



<n<! 



V2j 



\n — n' | 



t \n—n 



~2 



exp 



T 

(25) 

where the Associated Laguerre Polynomials L l n appear, 
with n< and n> denoting the smaller and larger, respec- 
tively, of the two indices n and n'. 



B. Solutions, difficulties, and limitations 



10 



200 
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10 
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10" 



Well 




7V=100 
N=200 
\ N=500 
• exp(7iw/2) 



100 200 300 400 500 



n 



From the explicit forms of the shift matrices in the 
two cases described in Eqs. (|2"5)) and (|2"5)) , it is clear that 
the shift matrix in Eq. (|18p . when inverted, has highly 
sin gular ele ments for N — > oo. To be precise, x n ~ 
n\J m\/m,2 for a square well (and yt n ~ y/2nmi/rri2 for 
oscillator) which implies that the elements in Eq. (|19l) of 
the shift-matrix have exponentially different scales. 

This difficulty of matrix inversion can be handled by 
performing a linear transformation from the set of basis 
states ip n (x) to a different set. Transformation to config- 
uration localized state is discussed in Ref. [18j. In our 
studies we used the singular value decomposition which 
is also effective. As was observed in Refs. (l9l . |25T - I28l | . 
there are complications in numerical convergence with 
an increased number N of included channels. The core 
of the problem is that the amplitudes for real and vir- 
tual channels involve very different scales. We find that 
for the square well, for example, remote virtual channels 



scale approximately as R n = R n Q ~ exp 

lar scaling follows for the oscillator) . The behavior of the 
virtual coefficients is illustrated in Fig. [T] Here and in 
what follows, the second index in R nn ' and T nn /, which 
corresponds to the incident channel, is dropped when it 
is the ground state channel, i.e., when n' = 0. 

As the physics of interest is comprised of small contri- 
butions from exponentially large excitations, the problem 
has mathematical issues. This is apparent also from Eqs. 
(jT6|) and ((20)) , where one is attempting to make a series 
with exponentially divergent coefficients vanish. This 
condition fails to fulfill especially for large mass-ratios 
7711/7712, when the coordinate-range \x\ ~ A (where ip n 
are not zero) implies large exponential factors e'' Kn '*' tl . 
The physics behind this is that when the interacting par- 
ticle is stopped by the wall, the non-interacting compo- 
nent continues its motion until its entire kinetic energy 



(simi- 



Figure 1: (Color online) Well and wall. Absolute values of 
reflection amplitudes for virtual channels are plotted against 
channel number n. The infinite square well potential (intrin- 
sic) is used in this example, with mi = 7712. The incident 
beam is in the ground state with kinetic energy E — 0.1. 
Different curves correspond to different truncations N. The 
straight line shows the curve exp(-7rn/2). The actual results 
closely follow this line initially and then deviate due to trun- 
cation in the channel space. 



is converted into virtual intrinsic excitations of the con- 
fining potential which is necessary for the system to be 
reflected. The bigger the mass of the non-interacting par- 
ticle mi relative to m2, the more kinetic energy it has, 
and the more complicated do the virtual excitations be- 
come. Figure [5] shows how this issue effects calculated 
results. 

For "good" mass-ratios, which is roughly when 
mi/m2 < 2, reliable solutions can be obtained (l8l. [20| 
that agree with exact and stable solutions gotten through 
a different method (the VPM, see Sec. IIVI) . as shown in 
Sec. IV Bl For these satisfactory results, the Projection 
Method had to involve arbitrary-precision numerics en- 
suring that both the small and the large contributions 
are properly taken care of. The reflection probabilities in 
different open channels calculated for the square well and 
harmonic oscillator models with m\jmi = 1 are reliable. 
They are also identical to those obtained through the 
VPM, and are shown in Figs. [T4l and [T51 respectively, in 
Sec. IV Bl These two figures display cusps at thresholds, 
which is a consequence of unitarity fl7L |29||. In addition 
to that, there are weak oscillations which, as discussed in 
the same section, become more pronounced in the case 
of a more massive non-interacting particle. 
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U n n'{X) (see Ref. [30]) through the following equation, 
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Figure 2: (Color online) HO and wall. This plot demon- 
strates the failure of the Projection Method, which uses Eq. 
(|18[) . to solve the scattering problem where the system of 
two particles bound by a harmonic oscillator confinement col- 
lides with an infinite wall. The incident kinetic energy is 
E = 1 = 0.5?^, which means that the total energy is half 
way between those of the ground state and of the first excited 
state. Different curves are labeled on the graph with the corre- 
sponding mass-ratios, mi/7712 = 1, 2, 3, 5, and 10. Phase shift 
S is shown to vary with N, the number of channels included 
in the calculation. The horizontal grid lines along with the 
tic-marks on the right indicate the values of the phase shifts, 
as obtained in a convergent way with a different method, the 
VPM (see Sec. IIV[l . The curve for mi/rri2 = 3 on a large 
scale up to TV = 160 is shown in the inset. 



The solution, however, is still elusive. This is espe- 
cially visible for "bad" mass-ratios with large values of 
mi/rri2. The Projection Method results are shown in Fig. 
[5] As the mass of the non-interacting particle mi gets 
larger, the results become extremely unstable. While a 
satisfactory value may be obtained for some cases, the 
approach is still flawed since the inclusion of more chan- 
nels (which must be accompanied by increased numerical 
precision) does not necessarily improve the results and 
may eventually lead to increasing oscillatory instabilities. 
This is demonstrated in the inset where the curve (3) for 
7711/7712 = 3, which seems to converge initially, is contin- 
ued up to N=160 where its behavior becomes erratic. 

In order to solve this problem, one should depart from 
projection onto the basis states. In what follows this is 
achieved by introducing a Variable Reflection Amplitude 



(26) 



Run' = e l ^ +K ^ x [2i^Kj^U nn ,(X) - 6 nn , 
Then Eq. (fTo) takes the form 

y^Unn' (HlX)ll> n (x) = 0, 



and can be efficiently solved by selecting a discrete set of 
N coordinate locations. This is the essence of a different 
approach discussed next. 



IV. THE VARIABLE PHASE METHOD (VPM) 

It follows from the discussion in the previous section 
that the approach based on projection onto the intrinsic 
basis is unpredictable in its ability to handle the prob- 
lem. As an alternative, the time-dependent methods have 
previously been used to treat similar problems [16| . Here 
we discuss the Variable Phase Method (VPM), which is 
a well established technique for treating multi-channel 
tunneling and scattering. It dates back to works pre- 
sented in Refs. [3ll437| . Exhaustive treatises on the sub- 
ject are found in books by Razavy [3(3], Babikov [38j|, 
and Calogero [3^]. Solving differential equations for the 
phases of the stationary-state wave functions, as func- 
tions of coordinate, is central to the VPM approach. 
These phases at asymptotic distances (from the scat- 
tering potential) make up the S-matrix of the problem. 
Equations for such quantities can be found by considering 
the phase shifts corresponding to the scattering potential 
being truncated at some coordinate locations. Alterna- 
tively, Green's function approach can be used. Tech- 
niques of this sort are widely used in reaction physics 
with atoms, molecules, and nuclei, and in relativistic 
scattering. Recently, there have been interests centered 
around multi-channel tunneling and scattering problems 
[3.l5,l3,l30l l40l - l4^ . Our problem is unlike those ordinarily 
encountered because its solution depends on proper treat- 
ment of the multi-channel virtual dynamics. Thus, here 
we enhance the VPM by applying it to virtual channels, 
which is mentioned in Ref. |38| as a possibility. Some 
later steps in this direction have been taken in Ref. [43| 
with off-shell amplitudes in the context of a three-body 
problem. 



A. Formulation of the VPM 

Let us first introduce the VPM briefly. We would like 
to emphasize that though we limit our discussion to one- 
dimensional scattering for simplicity, the approach is a 
general one. This method is also known as the Variable 
Reflection Amplitude Method [H, [H 0] . By using fac- 
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torization of the form 



for the wave function, Schrodinger's equation 
H$(X,x) = Et<&{X,x) with the Hamiltonian from 
([3|)-(j4l can be transformed into a coupled-channel equa- 
tion for the center-of-mass wave- functions ^ n (X) for 
channels n (subject to appropriate boundary conditions) 



d 2 



dX 2 



Kt 



,,{X)-^2v m >(X)^ n ,(X) = 0, (27) 



where the folded potentials are 
2M f°° 

V nn ,(X) = —j 1>*(x)V(X,x)4> n (x)dx, (28) 

and K n is defined in Eq. ([Tjl. 

The reflection and transmission amplitudes are spec- 
ified in reference to the potential-free solutions of 
Schrodinger's equation. These solutions are defined in 
terms of diagonal matrices as 



+iK n X 



y/-2iK n 



Inn' i 



(29) 



where the ± sign corresponds to a wave moving in the 
right/left direction. These solutions are normalized to 
unit current with the Wronskian set to unity, 



-E-(X) d ^l = l. 
y ' dX 



V ' dX 

The functions defined in (|29[) can be used for both open 
and closed channels provided that, as mentioned earlier, 
the principal branch of the square-root is selected for an 
imaginary K n . 

While there are variations of the VPM technique (30j . 
we demonstrate here the approach that explicitly empha- 
sizes the decoupling of the reflection and transmission 
coefficients and the different roles thereof [38]. 

It is convenient to apply the VPM by considering an 
auxiliary set of free-space wave functions 

V(X,X') = [E+(X) + E-(X) R(X')]T(X') (30) 

with coefficients R(X') and T(X') defined from the solu- 
tion Vt'(X) of Schrodinger's Eq. (f27"l) . using the Cauchy 
boundary condition at some point X', so that at X = X 1 



V(X,X') = V(X) and 



(31) 

It is convenient to interpret ^(X, X') for X < X' as 
the wave function corresponding to a potential truncated 
from the left, V nn >(X,X') = V nn >{X)6(X - X'), where 
9(X) is the Heaviside step function. So, V nn '(X,X') = 




Figure 3: An incoming wave of amplitude 1 is traveling from 
the left toward the potential barrier. It is reflected with 
amplitude R and transmitted with amplitude T. For any 
arbitrary point X' the barrier is thought of as a combina- 
tion of two parts: the unshaded part to the left of X' , and 
the shaded part V(X, X') to the right. At X' , in accor- 
dance with Eqs. (|30[> an d (|31|) . the incoming and outgoing 
components are identified as having amplitudes T(X') and 
R(X')T(X') respectively. Here T(X') represents the over- 
all amplitude of the wave function at X' which has been 
modified, relative to the incoming beam, due to the passage 
through the unshaded part of the barrier. In the context 
of the shaded part only, T(X') represents an incident beam 
normalization; thus R(X') is interpreted as the amplitude 
of reflection from the shaded part. Due to this normaliza- 
tion, the transmission amplitude T(X') through the shaded 
part is given by the final amplitude T (transmission through 
the full potential) normalized relative to the incident ampli- 
tude T(X'). Thus, T(X')T(X') = T. The incoming beam 
has a unit amplitude, and T(— oo) = 1, T(— oo) = T. It is 
obvious that R(—oo) = R which is the amplitude of reflec- 
tion from the full potential. On the right of the potential 
R(oo) = 0, T(oo) = T, and T(oo) = 1. 



at X < X', and the wave function is given by (f3"0")) with 
the boundary condition (|31[) . This interpretation is illus- 
trated in Fig. [3] Now, R(X') in (|30|) is a matrix in the 
channel space. With the help of Fig. [3l it can be iden- 
tified as the reflection amplitude for a wave scattering 
from the potential truncated from left (the shaded part 
in the figure). The vector T(X') in the channel space is 
the amplitude of the wave function ^(X, X') at X = X 1 . 
It can be normalized in different ways. In the litera- 
ture the wave function 1 F and its amplitude T are com- 
monly viewed as collections of independent column vec- 
tors corresponding to different independent initial con- 
ditions (initial channels). It is clear from Fig. [3] that 
transmission through the potential V(X, X') is inversely 
proportional to the amplitude T(X'). With relatively 
straight-forward derivations that involve substitution of 
the wave function in Schrodinger's Eq. (|2T|) by (f3"U)) (see 
also Ref. [38]), one can show that the matrix R(X) is 
subject to the differential equation 



dR(X) 
dX 



= [(5+ + R(X) 5-)] V [S+ + S- R(X)] . (32) 



The X-dependency of the S's and the folded potential V 
is suppressed in this differential equation and all others 
to follow unless there is an ambiguity. 
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The equation for the vector T(X) is linear, 



dT(X) 
dX 



= -E~V [E+ + E-R(X)] T{X), 



(33) 



which reflects linearity of quantum mechanics. To be 
more specific, the linearity shows that the column- vectors 
corresponding to different initial channels are indepen- 
dent of each other, and the amplitudes allow for arbitrary 
normalizations. 

The reaction physics of interest, in agreement with the 
discussion in Sec. Ill B[ is given by R alone; and, the 
physical properties are independent of normalization, re- 
sulting in the decoupled equation (j3"!S)) for R(X). 

As follows from the boundary condition on the wave 
function 9(X) or from the interpretation of R(X') as a 
reflection amplitude, R(X') is subject to the boundary 
condition 

-fiw(oo)=0, that leads to i?, m <(— oo) = R nn > ■ (34) 

The transmission amplitude is determined by T(X'). 
One can treat T(X') as a matrix of the column-vectors 
described previously. Assuming the incident beam to be 
in channel n' and normalizing it to unity (which is the 
most common and natural way), one would have 

T nn /(-oo) = S nn >, andthus, T nn > (oo) = T nn > . (35) 

The elements of the reflection matrix R(X) are sufficient 
to determine all observable probabilities. One, however, 
may still want to obtain the transmission amplitudes, 
which can be done by integrating Eq. ()33|) separately 
using previously determined values of R(X). Due to the 
different boundary conditions for R and T (see Eqs. 
and ()35[) ). this approach is computationally inconvenient. 
However, given that T ml '(oo) = T nn i, it can be inter- 
preted as a final state normalization to a yet-unknown 
value T nn i , and this inconvenience can be avoided. If one 
defines a matrix T(X) so that T(X)T(X) = T. Then 
T(X) coincides exactly with the matrix of transmission 
amplitudes through the truncated potential V(X, X') 
(see Fig. [3]). Since -^{T(X)T(X)) = 0, Eq. ^ in 
terms of T(A) is 



dT(X) 
dX 



= T(X)E~V [E++E-R(X)] 



It is to be used with the boundary condition 



T nn >(oo) = S n 



and then T nn >(— oo) = T n 



(36) 



(37) 



This approach is equivalent to the one discussed in Ref. 

In summary, the most celebrated advantages of the 
VPM are its physical transparency, generality, and the 
simplicity in its application. Phase equations can indeed 
be constructed for most quantum mechanical problems. 
In particular, with appropriate substitutions for the E 
functions, the approach can be immediately used in three 



dimensional problems with radial variables. The VPM is 
technically simple; the entire multi-channel problem is re- 
duced to a relatively straight-forward integration of Ric- 
cati equation (|3"!2)) from right to left with a zero starting- 
value (for R(X)) as a boundary condition. 



B. Virtual channels in the VPM 

In this work we find yet another value of the VPM in its 
effectiveness in treating virtual channels and in general 
complex- momentum (i.e., off-shell) applications. Some 
suggestions, in this direction, have been made in Refs. 
[38j, |43|. However, there are a few important things to 
note: 

First, the formalism remains valid in the complex mo- 
mentum plane assuming that for virtual channels the 
principal branch of the square root is selected in Eq. (JT]) . 

Secondly, separation of the amplitude (given by T(X)) 
from the physically relevant phase difference between in- 
coming and outgoing components (given by R(X)) is im- 
portant. Due to this separation, the principal equation 
(|3"!?|) is solved with a zero- value boundary condition (|34l) , 
and without any concern about exponentially falling or 
rising components of the wave function outside the po- 
tential. Normalization is provided by a set of decou- 
pled equations for each selected initial condition. There- 
fore, while solving scattering problems, one could con- 
sider only those columns T(X) that correspond to open 
channels of interest. 

Closed channels can be studied, if desired, with an ini- 
tial wave function exponentially rising toward the poten- 
tial. Bound states can also be explored in this way |39j, 
but we do not study these questions. 

Normalizing in a way to have the closed channels set 
to zero still deserves some attention. As further demon- 
strated in Sec. IVAl| T(X) and_T(V) for closed channels 
are set to zero differently. For T(X) one assumes the ini- 
tial beam normalization of zero for any closed channel, 
i.e., 



T nn '{— 00 ) =0 if n' is closed. 



(38) 



Thus after scattering one has waves in closed channels 
with amplitudes exponentially decaying to zero away 
from the potential. 

Contrary to that, eq. (|3"T)) is best thought of as a final 
state normalization, thus 



if n is closed. 



T„ n /(oo) = 



Finally, the original Eqs. ([3^)1 and (|3"3"|) used for open 
channels, are valid also for closed channels, but have is- 
sues with numerical stability for virtual excitations. The 
exponential divergence of the functions E^ n (X) with n 
(or, |-ftr ra |), especially at large \X\, makes it difficult to 
handle long-ranged potentials. We define 



U(X) 



E-R(X)E-, 
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which agrees with Eq. (|26l) . so that Eq. written in 

terms of the variables U(X) instead of R(X), reads 



dX 



■i{K n + K n ,)U nn , - UniVw U Vn , , (39) 



where U and V depend on X. It is noteworthy that in 
this form the equations no longer contain any exponential 
factors. A similar substitution can be done for T(X) or 

The amplitudes U for open channels oscillate outside 
the potential where V = 0, which is not the most desir- 
able boundary condition one would want to deal with. 
However, this is a minor inconvenience compared to the 
benefit of the exponential drop of U(X) for virtual chan- 
nels with distance from the potential barrier. This is 
particularly important because in the problems that we 
discuss, the reaction processes contain only a few open 
channels but are determined by numerous closed chan- 
nels. 



V. APPLICATIONS OF THE VPM 

A. A 5-barrier 

We proceed by considering a ^-barrier as the scatter- 
ing potential. Here, a bound system of two particles is 
incident on a potential 



V(xi, x 2 ) = 



AM 



8{x 2 ), 



(40) 



where, again, only the second particle interacts with the 
potential; A is a length parameter, characterizing the 
strength of the barrier. 

Interactions of composite objects, such as diatomic 
molecules, with a <5-barrier have been discussed before 
[30L |40L l4l| , but usually without any involvement of the 
virtual channels and in situations where the potential- 
barrier acts on both the particles. 

Any potential can be considered with the VPM ap- 
proach in principle; however, the short-ranged (^-potential 
provides a good way of exploring the generic features of 
scattering without putting efforts into computing folded- 
potentials. The folded potential (f28|) for a 5-barrier (140)) 
takes an analytic factorized form: 



Vnn'(X) = — C ( — 



Vv 



(41) 



where fii is the mass-ratio defined in (fl5j) . 

In the limit of A — > the (5-barrier turns into an im- 
penetrable wall, thus allowing us to complete the study 
in Sec. IIII1 which is done in Sec. IV Bl 

Introduction of the short-ranged 5-barrier adds just 
one additional length scale A to the parameters used to 
describe the problem. There are thus three length scales 



in the problem, namely, the intrinsic scale A, the inci- 
dent beam wavelength ~ 1/K, and the potential scat- 
tering length A. The corresponding energy scales are the 
intrinsic energy scale e, the incident beam kinetic energy 
E, and the energy scale associated with the (^-potential, 
defined by 



E s 



2MA 2 ' 



(42) 



The mass-ratio /x defined in Eq. ([15)) connects the 
length scales A and 1/K at similar energies; the precise 
relation is XK = WE/Jj,. 

The non-composite limit of the process is reached ei- 
ther if all the mass is concentrated in the interacting 
particle leaving fi± = and therefore A — > 00, or if the 
intrinsic states have infinitely high energy e — > 00 and 
thus A = 0. This yields a textbook problem of scatter- 
ing off a (5-barrier, where the transmission and reflection 
amplitudes are, respectively, 



T 



iKA 
iKA - 1 ' 



and R 



1 



iKA - 1 



(43) 



For a non-composite projectile the <5-potential allows for 
scattering only in the symmetric channel, since 



S + = 



1 + iKA 
I -iKA 1 



and S~ = 1. 



Transmission and reflection probabilities from a ^-barrier 
are determined solely by the energy ratio (KA) 2 = E/Eg. 
Thus, the sign of the coupling A, i.e., whether it is a well 
or a barrier, does not matter, 



\T 2 \ = 



E/E s 
E/E s + 1 ' 



|i? 2 | = 1- |T 2 | 



(44) 



1. Spatial dynamics of the reflection and transmission 
amplitudes 

The spatial dynamics of the reflection and transmission 
amplitudes is shown in Figs. f5]and[SJ As an example we 
take the "well" confinement (see Sec. IIII A 1[) with equal 
particle masses, fJ-i = fJ-2 = 1/2, and a 5-barrier with 
strength E$ = 1 in units of intrinsic excitations (|17|) . 
The kinetic energy of the beam is E = 4 in the same 
units, which means that there are two open channels. 

In each of these figures, the lower panel (c) shows 
folded potentials, as follow from Eqs. (|41[) and (|22l) . 
Thanks to the simplicity of 5-barrier; the folded poten- 
tials have obvious forms showing the structures of the 
wave functions for the intrinsic square well confinement. 
Naturally, V nn '(X) — outside the well, or, in other 
words, if \x\ > tt/2 (that is, \X\ > 7r/4) where x (as well 
as X) is expressed in units of A. 

As explained through Fig. [21 the dynamic transmission 
and reflection amplitudes, T(X) and R{X). in the VPM 
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-2 - \ / - 

-4 [ , , : 

-71/4 - 0.4 0.4 jc/4 
X 



Figure 4: (Color online) Well and 5-barrier. Barrier 
strength is E$ = 1. With incident beam kinetic energy 
E = 4, the only open channels are those labeled by or 1. 
The transmission and reflection probabilities, \T nn > (X)\ 2 and 
\R„„i (X)\ 2 , and the folded potentials V nn i(X) are shown in 
panels (a), (b) and (c) for the open channels. X is expressed 
in units of A. 

correspond to the potential truncated from the left of X . 
Therefore, both transmission and reflection probabilities 
\T nn i(X)\ 2 and \R nn t(X)\ 2 shown in panels (a) and (b), 
respectively, are evolved from right to left following Eqs. 
([53)1 and ![57]), For both Figs. Hand [5] we utilize the final 
state normalization (pTT|) of T(X) instead of T(X), due 
to the convenience in application and interpretation of T 
as the transmission amplitude. 

Figure [4] shows the dynamics in the open channels. 
Since the intrinsic potential is of a finite width, the fi- 
nal values of the reflection and transmission coefficients 
are reached at X = — 7r/4 which means inclusion of the 
full potential. And therefore, the values of \R nn i{X)\ 2 
and \T nn '(X)\ 2 at X = ±7r/4 are the asymptotic values 
thereof, commensurate with and (|57|) . The proba- 
bility is conserved at all values of X: 

\T nn '(X)\ 2 + \R nn ,(X)\ 2 = 1, 

nGopcn 

and R nn i(X) = R n > n (X) due to time reversal symmetry. 
However, T nn i(X) ^ T n i n (X) due to the asymmetry in 




-7C/2 -71/4 7C/4 

X 



Figure 5: (Color online) Well and 5-barrier. This graph 
refers to the same projectile as described in Fig. 2] with E = 
4. The reflection and transmission cross sections, |7? nn /(X)| 2 
and |T nn /(X)| 2 , the folded potentials V nn i{X) are shown for 
the out-going channel, n — 0, for a projectile in three virtual 
incoming channels n' = 2, 3, 5. X is expressed in units of A. 



the truncated potential. Symmetry is recovered in the 
final T since the full symmetric potential is covered at 
X = -7r/4, where \T nn ,\ 2 = |T„,„| 2 , as seen in Fig. gTa). 

A different picture emerges with the virtual channels. 
The linearity and independence of initial conditions, dis- 
cussed earlier, is important since it makes the normaliza- 
tion of virtual channels irrelevant for the S-matrix and 
other asymptotic reaction observables. Having said that, 
one can assume that T nn i{— oo) = if the initial chan- 
nel n' is closed (Eq. l[55|)). and thus T nn ,(X) = for 
any X due to linearity f|33[) . However, an incident beam 
in some open channel n' generates virtual excitations n 
that exist outside the potential. Therefore, T nn /{X) is 
an exponentially decaying non-zero function beyond the 
range of the potential when initial channel n' is open but 
the final channel n is closed. 

A totally different situation arises with the final state 
normalization (I37p . i.e., with T nn i (oo) = if n corre- 
sponds to a closed channel. Thus assuming that all vir- 
tual channels are normalized to zero to the right of the 
barrier generates disturbances of virtual channels in front 
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of the potential barrier, so that T nn i (X) is not zero there 
when n' is closed. This is seen in Fig. [3] where such 
quantities T nn t(X) (or, their norms) decay exponentially 
to the left of the barrier, but are not zero at X = —ir/4. 

The interpretation of R nn r(X) is quite different. It 
sets relations between the different components (phase 
shifts) of the wave function, and is non-zero for all real 
and virtual initial and final channels. Nevertheless, the 
behavior is similar (see Fig. [3]). 



2. Results 

Let us now discuss some final results for the scattering 
and tunneling of the deuteron-like system. Here we con- 
tinue to consider the infinite square well ("well") and the 
harmonic oscillator ("HO") models (see Sec. MI A II and 
IIII A 2|) . that do not allow for breakup. A model with a 
continuum of intrinsic states, which allows for breakup, 
is described in Sec. I VIII 

Figures [5] and [7] for the "well," and Figs. [5] and [5] for 
the "HO," depict the probabilities \T n \ 2 of transmission 
and | R n | 2 of reflection, respectively, in the first few lowest 
channels as functions of incident beam kinetic energy E. 
Each figure contains vertical grid lines indicating the lo- 
cations of channel thresholds. In the case of the "well," 
new channels open up at kinetic energies E — n 2 — 1 
where n is a positive integer; for the "HO," these occur at 
integral multiples of hui. In each case the incident beam 
is in the ground state channel, thus the corresponding 
subscript is suppressed, and n refers to the final chan- 
nel. The strength of the <5-barrier is set, via Eq. (|42|) . to 
E$ = 1. The redistribution of probabilities at the thresh- 
old energies, required by unitarity, leads to cusps in the 
cross sections (l7l . [29| . These discontinuities are common 
in all the figures for both the models. 

All figures contain three curves with fii = 
0.3, 0.5 and 0.7, and thus illustrates the mass-ratio de- 
pendence. In addition to that are shown the ground 
state to ground state transmission probabilities in the 
top panels of Figs. |B] and [5] in the non-composite limit 
/.ti — > where an analytic answer follows from Eq. (1431) . 
Indeed, when the non-interacting particle- 1 is very light 
compared to the interacting particle-2, i.e., when /ii — > 0, 
then particle-2 carries almost all the momentum, and the 
presence of particle-1 hardly matters. In this limit the 
behavior of the projectile approaches that of a single non- 
composite particle. 

The development of the resonant behavior, as [i\ in- 
creases, is easy to follow in these plots. For larger values 
of fii the curves exhibit prominent peaks and dips that 
are not associated with cusps at thresholds. Classically 
this can be viewed as a process in which the light interact- 
ing particle is stopped by the potential, while the larger 
mass fi\ keeps moving forward without any impediment, 
until most of its kinetic energy is transferred into poten- 
tial energy of the intrinsic interaction, and then it either 
turns back or pulls the smaller interacting mass through 
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Figure 6: (Color online) Well and 5-barrier. Barrier 
strength Eg = 1. Probabilities of transmission to different 
final channels n = 0, 1, 2 and 3 are shown as functions of in- 
cident beam kinetic energy E. The incident beam is in the 
ground state channel. The different curves on each panel cor- 
respond to different values of the mass-ratio fii , as labeled. 
The non-composite limit is shown in the top panel with the 
thick solid line. 
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Figure 7: (Color online) Well and 5-barrier. Same as Fig 
|6l showing reflection probabilities. 
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Figure 8: (Color online) HO and 5-barrier. Barrier strength 
Eg = 1 = 0.5 hoj. Probabilities of transmission to different 
final channels n = 0, 1, 2 and 3 are shown as functions of 
incident beam kinetic energy E. The incident beam is in 
the ground state channel. The different curves on each panel 
correspond to different values of the mass-ratio as labeled. 
The non-composite limit is shown in the top panel with the 
thick solid line. 
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Figure 9: (Color online) HO and 5-barrier. Same as Fig. 
[Sj showing reflection probabilities. 



Figure 10: (Color online) Well and 5-barrier. The system 
is same as that in Fig. [6] showing transmission probability 
|Tb| 2 for energies below the first threshold in panel panel (a), 
and phase shifts 8 for the same energy region in panels (b) 
and (c). 



the barrier. Hence, the larger the non- interacting parti- 
cle's mass, the more complex and chaotic the process. 

It is intuitive to suggest that the highly virtual chan- 
nels have little or no effect on observables at low-energies. 
It is proved otherwise in our studies. In Figs. HUl and [Til 
we focus on the low energy region, below the first thresh- 
old, for the "well" and the "HO" models, respectively. 
Here we use the same parameters as in Figs. and 
present our results for the non-composite limit fix = as 
well as for projectiles with /ii = 0.3, 0.5, and 0.7. Along 
with the transmission probability, we show results for the 
two phase shifts © that are defined up to the first thresh- 
old (shown by the vertical grid line). We conclude that 
the compositeness, and the composition of the projectile 
given by the mass-ratio of the components, are conse- 
quential factors that determine the observables. In these 
models, as well as in the ones with breakup (discussed in 
Sec. IVIIj) . we find a systematic enhancement of tunneling 
probability, with increasing mass of the non-interacting 
component, in a broad region of energy near the first 
threshold. This enhancement was earlier discussed in 
Ref. [20]. This is supported by a recent experimental 
at GANIL by Lemasson and others [Tl| that shows en- 
hancement in tunneling of heavy He isotopes, where ad- 
ditional spectator-neutrons contribute to the mass of the 
non- interacting component, while the alpha core inter- 
acts with the Coulomb barrier. 
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Figure 11: (Color online) HO and 5-barrier. The system 
is same as that in Fig. [8] showing transmission probability 
|T | 2 for energies below the first threshold in panel panel (a), 
and phase shifts 5 for the same energy region in panels (b) 
and (c). 



3. An attractive 8 -well 



Figure 12: (Color online) Well and attractive 5-well. The 

probability of transmission from ground state to ground state 
is shown as a function of incident kinetic energy. The upper 
and lower panels correspond to the barrier-strengths Eg = 1 
and 5 (i.e., 0.5 hw and 2.5 fkS), respectively. In both cases 
three different mass-ratios /Lti = 0.3, 0.5, and 0.7 are consid- 
ered, along with the analytic limit of a non-composite projec- 
tile (labeled by pi =0). 



In addition to the J-barrier discussed so far, we ex- 
plored scattering that involves an attractive 5- well. We 
stress again that for a non-composite projectile the sign 
of the interaction does not effect the observed reflection 
and transmission probabilities (see Eq. (|44| ). This is not 
true for a composite projectile. This topic has been ex- 
tensively explored in the literature and is often referred 
to as the Barkas Effect. In Ref. [45j, one can find more 
references that are relevant, and a model that is simi- 
lar in spirit and discusses the Coulomb excitation of a 
harmonic oscillator. 

Our results for the transmission probability in a scat- 
tering that involves a (5-well are shown in Fig. [121 We 
present results for the case of a harmonic-oscillator con- 
finement only; the results for the square well are similar. 
In all cases, even at relatively small masses of the non- 
interacting component, the scattering process is highly 
resonant. The interacting particle- 2 and the barrier form 
a bound state at an energy —^Eg, which is only a virtual 
binding in the three-body problem. However, the system 
in an excited state n with intrinsic energy e n can be tem- 
porarily bound as a whole, thus leading to a resonance at 
Et = £ n — fJ-2E$. We find that this crude interpretation 
unravels some of the complex resonant patterns seen in 
Fig. [T2] The resonances indeed periodically follow the 
channel thresholds, and they are close to the thresholds 



for small Eg = 1 (upper panel) and are further away for 
the larger Eg — 5 (lower panel). 



B. An infinite wall 

In this section we would like to return to the wall prob- 
lem which, as already shown in Sec. IIII[ is an extraor- 
dinarily illustrative example. This model emerges in the 
limit of a very strong ^-potential (Sec. IV Al) . i.e., with 
0. 

We study the convergence of the VPM method sepa- 
rately in Sec. IVH nevertheless, here we present Fig. Q2] 
which, in contrast to Fig. [21 shows that the VPM method 
is not prone to the convergence issues. Even for the large 
mass-ratio m\jmi — 5 the VPM produces a perfectly 
smooth curve that converges to a final 6 — —77°, which 
is not the case with the Projection Method. 

Figures [Til and [T51 show the reflection probabilities of a 
composite projectile in the ground state channel scat- 
tered from a wall. They are similar to the previous 
results for scattering that involves a 5-barrier (see Sec. 
IV A2[) ; cusps at thresholds and some resonant behavior 
are among the typical features. 

Scattering below the first threshold is characterized by 
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Figure 13: (Color online) HO and wall. This figure refers to 
the scattering problem described in Fig. [21 where the system 
of two particles bound by "HO" confinement collides with an 
infinite wall. The incident kinetic energy is E — 1 = 0.5 huj, 
and the mass-ratio is mi/m2 = 5. The phase S, as calculated 
using the VPM, is plotted in the solid line against the number 
of included channels N. The horizontal grid line indicates the 
final value of the phase shift to which it is found to converge 
smoothly with increasing N. The dashed curve shows results 
obtained through the Projection Method, which is unstable. 



a single phase shift 8, which is plotted in Fig. [16] as a 
function of incident beam kinetic energy, in the case of 
the harmonic oscillator. The different curves correspond 
to different mass-ratios. 

The limit of very low energies is particularly interest- 
ing. The formal effective range expansion (38l. l39l l46l l47l| . 
in the context of the VPM, has been applied extensively 
to problems of nucleon, molecular, and atomic scatter- 
ing. As K — > 0, the S-matrix, S — e 2lS , is characterized 
by a phase 8 — ~Ka where a is the scattering length. 
This length a depends only on the mass-ratio mxjmi 
and represents the distance of the turning point from the 
reflecting wall. A scattering length a > implies that the 
system is reflected at a distance a prior to reaching the 
wall. Figure 1X71 shows a in units of A, as a function of \x\. 
The limit \i\ — > corresponds to a non-composite case 
where the scattering length is zero. It is interesting to 
note that, while the intrinsic wave function of an infinite 
square well confinement has a finite width ~ 7rA, the scat- 
tering length can easily exceed this range. Thus, a clas- 
sically impossible situation occurs in which a finite-size 
system reflects from a wall before it actually approaches 
it within the contact distance. In the limit of fii — >1 the 
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Figure 14: Well and wall. This is the same as Fig. [7J but 

for a "wall," not for a 5-barrier, and only for the case of equal 
masses, mi = m2 (i.e., = 0.5). 
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Figure 15: HO and wall. This is the same as Fig. [9] but 
for a "wall," not for a 5-barrier, and only for the case of equal 
masses, mi = (i.e., /ii = 0.5). 











HO 





1 


i 


i 

1 










1 1 










1 


1 


1 





15 




Figure 16: (Color online) HO and wall. Phase shift as 
a function of incident beam kinetic energy. The curves are 
labeled with different values of the mass-ratio m\jmi- 



scattering length a diverges. This is a strong divergence 
since it is relative to a divergent scale, A — > oo for any 
given energy because of the vanishing reduced mass. It is 
worth pointing out that the divergence of the scattering 
length due to intrinsic degrees of freedom coupling to the 
reaction dynamics is known as Feshbach Resonance. 

For both square well and oscillator models one can ex- 
amine the analytic results for a by considering a few vir- 
tual channels within the Projection Method. It becomes 
immediately clear that such an expansion is convergent 
only in the limit of [i\ — > 0. In this limit we obtain 
a/A « 0.56 /i^ 2 for the square well bound system, and 
a/A = for the oscillator-bound system. 

Figure [TH] shows the square of the amplitude of the 
wave function, £|T n „< (X)| 2 , for a projectile in the in- 

n 

coming channel n' = 0. This is interpreted as the density 
of probability for the center of mass of the projectile to 
be at a location X when it is reflected from an infinite 
wall. The four curves show a few of the most represen- 
tative situations; incident beam kinetic energies as low 
as E — 1.5 and as high as 30, and two different mass- 
ratios Hi = 0.5 and 0.9. All probabilities eventually die 
to zero beyond the wall located at X = 0. The first two 
curves represent cases where the energy of the projectile, 
E = 1.5, is halfway between the energies of the ground 
state and the first excited state. Hence only one open 
channel is present. For mass-ratio /ii = 0.5 the behav- 
ior is plain. However, when the non-interacting particle 
contains 90% of the total mass there is a peak of proba- 



Figure 17: (Color online) Well/HO and wall. For the 

"deuteron and Coulomb wall" model, the scattering length 
a (in units of intrinsic length A) is shown as a function of 
the mass-ratio /ii = mi/M for two different systems: those 
bound by the "well" and by the "HO" confinements. 

bility density in front of the wall. This is consistent with 
the enhanced scattering length (see Fig. \T7} and with 
its interpretation that this probability peak corresponds 
to a turning point where the system is stopped prior to 
reaching the wall. At higher beam energies the center 
of mass penetrates considerably through the wall (region 
X > 0). As expected, this penetration is deeper for a 
more massive non-interacting component; the peaks in 
the density inside the wall can also be attributed to the 
non-interacting particle being stopped via energy trans- 
fer to intrinsic excitations. 



VI. ROLE OF VIRTUAL CHANNELS AND 
CONVERGENCE 

While problems similar to those presented here have 
been extensively discussed in recent literature, for exam- 
ple Refs. 0-0, little attention has been paid 
to virtual channels. In fact, most of these works discuss 
tunneling of a diatomic molecule where both the atoms 
interact with the potential. In that case, virtual excita- 
tions are relatively less likely to take place, and hence the 
folded potential within open channels already provides a 
relatively good description of the process. Our selection 
of models on the other hand, where only one particle in- 
teracts with the scatterer, is dynamically different. It is 
the virtual channels that shape the non-interacting par- 
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Figure 18: (Color online) Well and wall. Probability den- 
sities as functions of location, for scattering of a square-well- 
bound system off an infinite wall at X — 0. 



tide's movement. Therefore, compared to the models 
discussed by other authors cited above, our models are 
in general more sensitive to virtual channels. In the most 
extreme case of reflection from an infinite wall, no mean- 
ingful description is possible at all without reference to 
the virtual channels. The folded potential for the ground 
state, depicted in Fig. U^c) with the solid black line, has 
a single hump, and therefore does not lead to any res- 
onant behavior in reactions at low energies, when only 
one channel is open. Hence it can be concluded that, 
the resonance-like increases or decreases in the transmis- 
sion and reflection probabilities shown, for example, in 
Figs. [5] and [51 at low energies and especially when the 
non-interaction particle is heavy, are exclusively due to 
virtual channels. 

A successful extension of the VPM so as to include 
virtual channels in the formalism, and the study of their 
role that we discuss in this section, are among the main 
achievements of this work. 

The importance both of compositeness and of virtual 
channels is illustrated in Fig. Q1J1 where we consider an 
oscillator scattering from a <5-barrier. The curves in three 
different styles and colors correspond to results from 
three different calculations: the solid red line represents 
the exact solution, i.e., the solution of scattering of a 
composite projectile obtained through a calculation that 
includes the virtual channels as well as the open channels; 
the dotted blue line represents scattering of a composite 
projectile solved through a folded potential but ignoring 
any virtual channel whatsoever; the solid black line repre- 
sents scattering of a non-composite projectile of the same 
mass as that of the composite projectile. For the region 
of energies shown in these graphs, there is only one open 
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Figure 19: (Color online) HO and 5-barrier. A two-body 
system with mi = mi, bound by an oscillator confinement, is 
scattered by an external delta potential with Eg = 1. Panels 
(a), (b), and (c) correspond to transmission probability, and 
two symmetric and antisymmetric phase shifts, respectively. 
Three curves shown compare full solution (labeled "full") with 
the approximate treatment that includes only the open chan- 
nels ("open"), and with results for non-composite projectile 
("non-composite"). We conclude that ignoring the composite 
nature of the system or neglecting the virtual channels results 
in neither the phase shifts nor the transmission probability 
correctly. 



channel; thus the asymptotic behavior of the wave func- 
tion is fully determined by the two phase shifts 6 which 
are shown in panels (b) and (c) as functions of incident 
beam kinetic energy. All three curves are different, indi- 
cating that neither non-compositeness nor treatment of 
open channels only can substitute for a full solution. To 
emphasize this, we show in panel (a) transmission proba- 
bility, an observable quantity, as a function of energy; for 
most of the energy region shown, the actual transmission 
probability appears to be higher than that of an equally 
massive non-composite particle. 

While the virtual channels cannot in general be ig- 
nored, the contribution of the highly excited states is ex- 
pected to diminish. Practical applications require some 
truncation in the channel space, too. In Fig. I2"07a) we 
demonstrate the rate of convergence of the transmission 
probability Pt — \Tq\ 2 by plotting it as a function of the 
number of included channels N. The curve is visually 
indistinguishable from the hyperbola 



P T (N) = P T 



Ns 
N 



(45) 
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Figure 20: (Color online) Well and <5-barrier. Convergence 
of the transmission probability Pt = | To | 2 to its limiting 
value, as a function of N. The masses are equal, mi = mi, 
and the kinetic energy is E = 4, so that only the first two 
channels are open. The strength of the S- barrier is Eg = 1. 
Plot (a) shows the actual behavior of the transmission prob- 
ability Pt(N) with N, and its asymptotic value Pt = 0.6708 
shown by the grid-line. In order to show agreement with Eq. 
dHJ, the value N [P T - P T (N)] is shown in plot (b). This 
quantity is well described by a constant N s — 0.157, shown 
with the red horizontal line. 



showing that the deviation of the probability Pt(N) from 
its limiting value Pt is inversely proportional to N. That 
is, N s is the rate of convergence. The lower plot, Fig. 
[207 b). shows agreement with Eq. (|45|) by comparing 
[Pt — Pt(N)] N with a constant N s . These results are 
for the square well bound system with \Kjf\ ~ N. With 
more precise consideration it is found that in general, 
the amplitudes converge as ~ 1/\Kn\. Figure demon- 
strates an excellent agreement with this rule using two 
different systems reflecting from an infinite wall. 

This power-law convergence is slow in contrast to an 
exponential convergence usually encountered for eigen- 
values and other structural observables as functions of 
truncation [24j . Here we repeat our recent conjecture 
[j"o| that this is an inherent property of reaction physics, 
where the kinetic energy operator plays a major role in 
the Hamiltonian. The mentioned operator discretized 
in coordinate space corresponds to a tri-diagonal matrix 
that meets a set of criteria for the power-law convergence 

In the course of our work we have vigorously tested 




0.01 



10 



100 



N 



Figure 21: Well/HO and wall. Convergence of the phase 
shift. The log- log scale acts to accentuate a good agreement 
with the ~ 1/jA'jvj convergence rule. The dashed straight 
lines represent these rules. For the square well intrinsic po- 
tential, K N ~ N and S(N) ~ 5 - N s /N with 5 = -23.05° 
and N a — 11.3. For the harmonic oscillator Km ~ y/N, so 
5(N) ~ 5 - y/WT/N with 8 = -22.98° and N s = 53.4. For 
both cases the incident beam energy is halfway between those 
of the ground state and the first excited state. 



the ~ 1 / 1 -ftTjy | convergence rule. While the rate of con- 
vergence, N s , depends strongly on the type of the system, 
we found no exception from the power-law convergence. 



VII. INTRINSIC POTENTIAL WITH A 
CONTINUUM; BREAKUP 

We would like to conclude our exploration with a some- 
what more realistic situation where the intrinsic potential 
allows for a breakup. We therefore consider a confining 
potential v(x) that has both bound state(s) and a con- 
tinuum. In our discussion below we study a particular 
confinement, namely, a finite square well ("finite well"): 



v(x) 



when |x| > A 
—v otherwise 



(46) 



This allows one to capture the generic features of the 
problem, while still having a small number of parameters 
and retaining the ability to have analytic solutions ([S]) for 
the intrinsic Hamiltonian. Realistic applications to three- 
dimensional problems with other potentials are outside 
the scope of this work. 

In what follows we again select the units of length A to 
represent the width of the intrinsic potential as defined in 
Eq. (|46l) . The bound state energies for the finite square 
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Figure 22: (Color online) Finite well and 5-barrier. Con- 
vergence of the transmission probability to its final value with 
truncation energy (in units of e) of virtual excitations. The 
well-depth is v — 1 and the incident kinetic energy is E — 0.2. 
Therefore, only the elastic channel is open. The figure con- 
tains four different curves, three of which are labeled with 
the widths of the quantization box, L = 15, 50, and 100 that 
were used. These curves agree well with the power-law con- 
vergence discussed in the previous section, the fit for which 
is shown with the last curve (green, dashed). The final value 
for transmission probability is 0.174. 



well potential are given by the transcendental equation 

±1/2 



tan \Je + v = ± 



v + e 



(47) 



where we remind the reader that v and e are expressed 
in units of e [see (|17[>]. The ± sign corresponds to the 
intrinsic parity V = ±1 of the state of interest. In addi- 
tion to these bound states there is a continuum of states 
with positive energies above the well. 

To model such a situation mathematically we discretize 
the spectrum using a quantization-box of width 2L, and 
hence the intrinsic wave functions ip n (x) are subject to a 
boundary condition ip n (±L) = 0. The choice of a large 
enough L can yield a spectrum that represents the con- 
tinuum obtained without the box. The large box allows 
for any finite potentials v(x) to be considered. 

In order to examine the appropriateness of the ap- 
proach and to address the potential concerns arising from 
the presence of the continuum and its truncation, we 
show in Fig. [55] the calculated transmission probabil- 
ity for such a system incident on a <5-barrier. A depth of 
v = 1 has been chosen for this example, thus there being 



a single bound state at energy Eq = —0.454, as follows 
from Eq. (g7|), with the RMS size 1.17 of the wave func- 
tion. The incident kinetic energy of the center-of-mass 
motion is assumed to be E — 0.2, which means that only 
the elastic channel is open. Though there is not enough 
beam energy for a breakup, the virtual channels are still 
important. In Fig. [55] we explore different box-widths, 
L = 15, 50, and 100 (in units of A), which shows that the 
results are independent of L, if it is large enough. (Note 
that with a large box-width the density of states in the 
continuum is high and therefore it is difficult to include 
high-energy channels.) Even for the smallest box L = 15 
the energy of the ground state differs from the exact an- 
swer only by less than 0.05%. The breakup threshold is 
at 0.495, which is slightly different from — Eq mainly be- 
cause the first excited state (continuum threshold) in the 
box does not exactly coincide with zero energy. These 
differences are minor and orders of magnitude smaller 
for L = 100. 

Now that the appropriateness and validity of the ap- 
proach is established, we present the transmission prob- 
ability and both symmetric and antisymmetric phase 
shifts S ± as functions of incident kinetic energy in Fig. 
[551 This figure concentrates on the energy region below 
the breakup threshold. This situation is important, since 
it is commonly encountered in practice. As it is clear from 
the graphs, the composite nature of the system and the 
virtual continuum are playing a crucial role in shaping the 
reaction process. We find that at very low energies trans- 
mission is inhibited for a composite particle. At higher 
energies close to the breakup threshold, transmission rate 
is always enhanced. Moreover, this rate increases for an 
increasingly heavy non-interacting particle. The role of 
the virtual channels appears to be universal for all mod- 
els that we investigated [see also Figs. [5] and [5], where 
|To| 2 increases sharply near the first threshold. The ex- 
periment in GANIL, as cited at the end of Sec. IV A 21 
proves tunneling enhancement for systems with breakup 
also. 

It is interesting to review the distribution of probability 
of breakup into the continuum, when energetically pos- 
sible. In the continuum we can still separate the center- 
of-mass and the relative kinetic energies. The relative 
kinetic energy is now given by the discretized states in 
the box. Fig. [541 shows the normalized probability distri- 
bution for breakup with different relative energies. The 
initial beam in this case corresponds to a projectile in 
the ground state with kinetic energy of E — 2. Thus, 
the total kinetic energy of fragments after the breakup is 
E + Eq — Ei ~ E + £o = 1-55 which is viewed as a sum of 
two parts: the center-of-mass kinetic energy h 2 K 2 /(2M) 
and the relative kinetic energy e. As seen from the plot, 
it is most likely to have the two fragments moving to- 
gether with very little relative energy (corresponding to 
the peaks on the left side of the plot), or, inversely, mov- 
ing apart in opposite directions with most of the energy 
concentrated in the relative motion (corresponding to the 
peaks on the right side). 
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Figure 23: (Color online) Finite well and 5-barrier. Trans- 
mission probability and phase shifts for the scattering of a 
two-particle system with intrinsic potential that is finite and 
allows for breakup. The potential is modeled by a square well 
of unit width and depth v = 1. We concentrate on the re- 
gion below the breakup threshold of kinetic energy at 0.495 
shown with the vertical grid-line. Four curves correspond to 
four different mass-ratios \x\ — 0.3, 0.5, and 0.7 and the non 
composite limit of for which Eq. (|43[) is plotted. The 5- 
barrier strength is assumed to be Eg = 1. The upper panel 
shows transmission probability; the curve is continued above 
the breakup threshold to show the cusp at the threshold. The 
two lower panels, similar to those in Figs. llll and [T0l show the 
phase shifts 8 + and S~ as defined in Eq. [5] Note that their 
meanings as phase shifts of the S-matrix is true only below the 
breakup threshold; above the breakup these are arguments of 
the corresponding reflection and transmission amplitudes. 

The choice of quantities plotted in Fig. [53] reflect our 
method, but is not very closely related to a potential ex- 
periment, where the momenta or velocities of both par- 
ticles could be measured. Therefore in Fig. [25] we show 
the distribution of probability to observe a certain com- 
bination of the particle- velocities vi and V2- Since the 
total kinetic energy after breakup is fixed, 

i (m 1 v? + m 2 v|) =E + e , (48) 

it is sufficient to use a single angle to parametrize the po- 
sition on the ellipse formed in the velocity (vi, V2) plane. 
In Fig. [23 we show the probability as a function of an- 
gle using a contour plot. Both Figures [24l and [25l have 
been smoothed, but preserve the general shapes which we 
believe to be good representations of the physics. Some 
of the features seen in Fig. [M] become more transpar- 



Figure 24: (Color online) Finite well and 5-barrier. Prob- 
ability to break up into a level e in the continuum is plotted 
as a function of e. The intrinsic potential is modeled by a 
square well of unit width and depth v = 1, as described in 
the text. The kinetic energy of the incident wave is E — 2. 
The system interacts with a <5-barrier of strength Eg = 1, 
and then either scatters elastically, or breaks up into its con- 
stituents with relative energy e. The three curves correspond 
to three different mass-ratios /ii = 0.3, 0.5, and 0.7. The cor- 
responding probabilities for elastic scattering are 0.64, 0.60, 
and 0.69 respectively, which complement the breakup proba- 
bilities shown here. 

ent in Fig. [25] For all mass-ratios we see that the 
probability peaks at about —45° which corresponds to 
non- interacting particle- 1 moving forward and interact- 
ing particle-2 being reflected back with velocities nearly 
equal in magnitude and opposite in direction. For equal 
masses the peak at low relative kinetic energies in Fig. EMI 
appears as two peaks in Fig. f2S]at about 45° and —135°. 
In both cases the particles move with similar velocities, 
forward for 45° and backward for —135°. The observed 
picture appears to be quite intuitive. 

VIII. SUMMARY AND CONCLUSIONS 

In this work we revisit some of the most intricate 
questions of reaction physics involving composite objects. 
The research presented here was inspired by highly un- 
predictable behavior of reaction observables including 
resonances and cusp-discontinuities, a very broad spec- 
trum of scales involved, and, at the same time, the ut- 
most importance of these processes in nuclear physics and 
other fields. 
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Figure 25: (Color online) Finite well and 5-barrier. The 

system and scattering conditions are the same as in Fig. 1241 
The system interacts with the 5-barrier, and then either scat- 
ters elastically, or breaks up into its constituents with mo- 
menta Vi and V2. Break-up probability, as a function of ve- 
locities vi and V2, is plotted as a contour graph. Due to Eq. 
(|48[> the points representing the possible velocity values lie on 
an ellipse in the (vi, V2) plane and are parametrized by a sin- 
gle angle. The contour plot shows the probability of breakup 
(radial coordinate) as a function of this angle. 



The topic of reactions involving composite objects is 
widely investigated, and benefits from many advanced 
methods and techniques. Our work got its thrust from a 
simple and well-defined problem of a deuteron-like system 
interacting in one dimension with an infinite Coulomb 
wall. A number of methods tried, in the past, for this 
problem have failed under certain circumstances, namely, 
where approximations or simplifications commonly used 
were not appropriate for this particular problem, or 
where the method did not produce convergent results, or 
where there were difficulties with numerical errors. Due 
to this delicate nature of the problem, we find exact so- 
lutions to all the examples considered in this work. This 
allows us to obtain comprehensive answers, and to be 
able to carry out comparisons with other methods and 
solutions. 

This goal requires us to have a precisely-defined Hamil- 
tonian with as few parameters as possible, and condi- 
tions that perhaps are more critical to reaction-structure 
interplay than those typically encountered in nature. 
In this presentation we restrict our discussion to mod- 
els. Nevertheless, our methods have broad applicability; 
most examples can be modified easily to represent re- 
alistic situations, and we continuously suggest cases in 
nature that are similar to what we discuss. We study, 
through this work, a two-particle system interacting in 
a one-dimensional scattering with a target that poses a 
^-potential or an infinite wall potential. It is always as- 
sumed that only one of the two components interacts 



with the target. The study includes models that do al- 
low the projectile to breakup, and models that do not. 
The dominant and non-perturbative role of the virtual 
channels that extend far in excitation energy is the main 
common theme of all the examples discussed here. 

We start by revisiting the "deuteron and Coulomb- 
wall" model which has been discussed for almost a decade 
with little outcome [18|, |20j, [25J . Unfortunately, this prob- 
lem is commonly dismissed either at the first glance when 
it seems uninteresting, or after some investigation when 
it seems unphysical, ill-defined, or unsolvable. We, on 
the other hand, find this model remarkable in its ability 
to demonstrate, in an extremely transparent manner, the 
dynamics driven by the virtual excitations. 

We review, and carefully apply, the technique of pro- 
jecting the reaction dynamics onto an intrinsic space and 
show that while satisfactory results are obtained in some 
limits, this formally exact approach does not yield con- 
vergent solutions in general. This is an important finding 
because this "Projection Method" is a prototype of sev- 
eral commonly used approaches in many-body problems 
that involve both structure and reactions [2ll |22||. 

As our main workhorse we utilize the Variable Phase 
Method (VPM) to address the coupled-channel problems 
of interest. While the method has been used by others 
before, we modify and extend it to treat highly remote 
virtual channels. We demonstrate that the VPM pro- 
duces reliable and convergent results. We investigate the 
contributions from remote virtual excitations, study con- 
vergence with the truncation size, and find the power-law 
convergence, which is in contrast to the exponential con- 
vergence seen in many-body structure problems [24j]. 

Within a given set of models, this work contains 
numerous examples, investigations and demonstrations. 
Cusps and discontinuities appear in observables as man- 
ifestations of conservation of probability and redistribu- 
tion of flux at the thresholds. Intrinsic structure gives 
rise to resonance-like behavior in tunneling probabilities; 
our models and recent experimental evidences indicate 
a generic enhancement in transmission probabilities due 
to virtual channels or a virtual continuum, whichever 
is the case. We explore and discuss the role of virtual 
excitations at very low energies, showing that even in 
those cases the scattering length is sensitive to the pro- 
jectile's structure. Due to the intrinsic structure and its 
coupling to reaction dynamics, scattering length can be- 
come infinite, the phenomenon being known as Feshbach 
resonance. We demonstrate how the intrinsic structure 
violates charge symmetry, which is called the Barkas ef- 
fect. The scattering of a non-composite projectile off a 
<5-barrier is the same for attractive and repulsive interac- 
tions. But, in case of a composite projectile, the corre- 
sponding three-body problem for an attractive potential 
is quite different from that for a repulsive barrier, and 
reveals numerous resonances, some of which can be un- 
derstood as bound states built upon individual intrinsic 
excitations involving two-body subsystems. 

Scattering and breakup dynamics influenced by a vir- 
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tual continuum are also investigated in this work. It is 
seen that the most probable breakups take place where 
either almost all the kinetic energy is relative, or almost 
all of it is in the center of mass. 
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